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Abstract 

> 

£\j . A formulation for performing quantum computing in a projected subspace 

o, 

is presented, based on the subdynamical kinetic equation (SKE) for an open 



X 



quantum system. The eigenvectors of the kinetic equation are shown to re- 
main invariant before and after interaction with the environment. However, 



the eigenvalues in the projected subspace exhibit a type of phase shift to the 

C3 , 

evolutionary states. This phase shift does not destroy the decoherence-free 

CT 1 ' 

(DF) property of the subspace because the associated fidelity is 1. This per- 
mits a universal formalism to be presented - the eigenprojectors of the free 
part of the Hamiltonian for the system and bath may be used to construct a 
DF projected subspace based on the SKE. To eliminate possible phase or uni- 
tary errors induced by the change in the eigenvalues, a cancellation technique 
is proposed, using the adjustment of the coupling time, and applied to a two 
qubit computing system. A general criteria for constructing a DF projected 
subspace from the SKE is discussed. Finally, a proposal for using triangula- 
tion to realize a decoherence-free subsystem based on SKE is presented. The 
concrete formulation for a two-qubit model is given exactly. Our approach is 
novel and general, and appears applicable to any type of decoherence. 
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I. INTRODUCTION 



Solid-state proposals for quantum computers have progressed markedly recently: these 
include superconducting junctions, quantum dots, electron spin resonance, nuclear spins of 
impurity atoms, and nuclear spins in a crystal lattice' 1-5 ^. The working conditions for most 
proposals require low operating temperatures (typically a few K). Ideally, if such system can 
be created as isolated from the environment, system evolution for a given quantum comput- 
ing process, may be described by a unitary operator with time reversal symmetry. Often this 
is not the case for real quantum computing systems because the interaction between such 
systems and their environment introduce decoherence, which destroys the superpositions of 
qubits that enable the quantum logical operations to be validated' 6 ' 7 ' . Such decoherence is 
a major obstacle in developing a practical quantum computer. Recent publications have 
formulated a theory for a decoherence-free (DF) subspace in which quantum computing is 
performed. The first formulation of DF conditions was performed by Paolo Zanardi and 
Mario Rasetti' 8 ' 9 '. This original formulation is quite general, which is not in the form of 
a master equation and does not invoke the Born-Markov approximation. Further work by 
Viola and Zanardi deal with DF subspaces by quantum control theory in a non-Markovian 
setting' 10 ' 11 ' 12 '. Then an important direction of the DF subspace theory was developed from 
a Master equation (such as the Lindblad equation) for the open system within the powerful 
semigroup approach' 13-17 '. A quantum jump description for the state of the atoms has also 
derived a consistent result with Master equation for dipole interaction atoms, but may give 
more insight into the time evolution of a single system' 18 ^. Experimentally, DF subspace 
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have recently been observed, which shows that such DF subspaces do indeed exist, allowing 
logical qubits to be encoded without decoherence^ 19 ' 20 ]. These provide the motivation for 
building appropriate open quantum computing systems suitable for (a) canceling the effects 
of intrinsic decoherence, and thus behaving as ideal quantum computing systems in the 
subspaces, or (6) construct DF subspace without invoking any approximations (such as the 
Born-Markov approximation) or restrictions on the type of decoherence (e.g., symmetric and 
collective decoherence). In this open system, self-adjoint operators and unitary evolution 
groups are not intrinsically required for quantum computation. Quantum computation can 
then be performed in a more general functional space including triplet structures: Rigged 
Hilbert Space (RHS), $ C H C $ x , or Rigged Liouville Space (RLS), $ ® $ GH®H C 
x $ ® $ x , where H is Hilbert space, $ is dense subspace of H and $ x is dual space of <3>, 
rather than just Hilbert or Liouville Space' 21-23 '. 

In this work, we present a formulation for performing quantum computation in a DF 
projected subspace based on the subdynamical equation (SKE). In section II, we briefly 
introduce a subdynamical formulation in the Hamiltonian representation. In section III, we 
propose a type of DF projected subspace. In sections IV and V, these concepts are applied 
to the example of a two qubit spin computing system plus Bosonic bath; a procedure to 
cancel the phase error induced by the change of eigenvalues is discussed. In section VI, 
a general approach for obtaining the necessary and sufficient condition for DF behavior is 
discussed. In section VII, a triangulating method is introduced for determining the DF 
projected subspace exactly. 

II. SUBDYNAMICS FORMULATION 

Consider a quantum system which is composed of iV quantum registers and interacts with 
a large thermal reservoir. We denote H$ (t), Hb and H int as the Hamiltonian of the system 
S, the Hamiltonian of the the thermal reservoir B, and the interaction between S and B, 
respectively. Then the total Hamiltonian of the system plus the reservoir can be expressed 
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as H s (t) ®Ib + Is® H b + XH int , with the corresponding Schrodinger equation as (t) = 
if (t) ip (t). If one chooses the time-independent eigenprojectors of H s (t) <g> I B + Is ® ife as 
P v = |<^,) and Qi, with Pj, + Qi, = 1 and the eigenprojectors of the total Hamiltonian as 
Yi v (t) to satisfy the Heisenberg equation, then, from the definition of U u (t) one can induce 
the definition of the creation correlation operator C v (t) (D v (t)) as Q u Ti u (t) = C v (t) (t) 
(IT, {t) Q v = U u (t) D u {t)), respectively. Note that the operator C u {t) creates the Q^-part 
of Ii u (t) from P v and the operator D v (t) destroys the Q^-part of Tl u (t) from P V) since 
C v (t) = Q V C V (t) P u and D v (t) = P V D V (t) Q v . This enables a projected kinetic equation to 
be constructed in the projected subspace by 

4 £ (*) * (*)) = °WE (*) ^ (*)) , (i) 

with an intermediate operator defined as 

6 (i) = (if 5 (t) + P B + ff* n ') + (iff (*) + iff + if*"*) C (i) (2) 
= O (f) + H?* + (iff (*) + iff + if*"*) C (t) , 

here JjS,B,mt ^ j^s,B,mt ^ arg ^ e diagonal ano i off-diagonal parts of the corresponding 
total Hamiltonian, respectively. 11^ (t) can be found from the subdynamics formulation as 
{P v + C v (t)) {Pv + D v (t) C v {t)Y x {P v + D v (t)) [21,22] . The creation (destruction) operator 
can be obtained from the basic operator equations in the subdynamics formulation. Indeed, 
from the definition of the eigenprojectors Ii v (t) we have 

{E u (t) - Q U H (t) Q v ) Q U U U (t) = Q u H int P u U u (t) , (3) 

and hence yielding the crucial relationship 

Q u U u {t) = Q u C{t)P u U u {t), (4) 

with definitions 



and 



C v (t) = Q v C(t)P v (6) 



" E v (t) - Q V H (t) Q v 
where E v (t) is z/th eigenvalue of the total Hamiltonian H (t). In the same way, 

D v (t) = P ^ntQ VEAt) _ Qi/H{t)Qv - (7) 

If one defines a projected wavefunction, ^ pro - ? (t), by ip (t) = (Pv + C u (t)) ^) V J°^ (t), then 

V 

J2 U Pu^-u (t) ip (t) = Yju i^J^ (t) = i ,proj (t), giving the formal evolution equation for a pro- 
jected wavefunction ip pro i (t) as 

^proj ^ = f e -iJ t0 dt'B(t')^ proj ^ ^ ^ 

= fe 1 ^ dt ' e{t,) ^ (Pu + D v (t ) a (to))" 1 (P, + A, (to)) V> (to) • 



III. DF PROJECT SUBSPACE 

Using the above subdynamics formulation, it is apparent that before the interaction 
occurs, the spectral decomposition for (H s (t) + H B ) is the same as that for the total inter- 
mediate operator (t), i.e., 

H s (t) + H B = Q(t) = J2 Pi (*) Pu, (9) 

V 

where E® (t) is the z/th eigenvalue of (H s (t) + H B ). But once the interaction with the 
environment has occurred, the eigenvalues change from E® (t) to E® (t) + AE l/ (t), where the 
definition 

AE v {t) = {(j> v \Hi ra C{t)\(f> v ), (10) 

is made. The corresponding eigenprojectors P v of (t) still remain invariant. The spectral 
decomposition of the intermediate operator (t) thus is different from (Hs (t) + H B ), 

© (t) = £ {K (t) + AE V (t)) P v . (11) 

V 
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The corresponding mixed-state fidelity in the projected subspace can be calculated from 



F(t) = Tr^pP r °i (t ) pP roj (t) 



(12) 



Tr. 



x , E (P pr ° J (*o))„„ ^) ff e"*/^ ^ E (*>))„„ PvrfjfW" 
E [(p" roj (to)) ^e-if^WW+^VW 'p^QW+^t'M' 



,i e [(^ (*»)) j^^;iw'**'* ft /rj*)«^* m 



= E V[(P proJ (to)) J = Tr {\jP pr ° 3 (to) P pr ° 3 (to)) = I- 

where (ff ro: > (to)) are the matrix elements of the density operator p pro - ? (to)- This exposes 
an exciting result: there is no decoherence in the projected subspace for states of the system 
since the eigenvectors remain invariant. However, a change in the eigenvalues introduces a 
phase shift in the evolution of the states, 

e -ifEi(r)dr ^proj ^\ Af ter the interaction e -/(^W+AB,W)^ ^ ^ 



From these findings there are two conclusions: 

(1) In general, for any system S + B, we can use the eigenprojectors of H s to construct 
a DF projected subspace in which the eigenprojectors remain invariant before and after the 
interaction between S and B, while the eigenvalues induce a phase shift in the eigenstates. 
The encoded states in the projected subspace are the projected states which are related to 
the original states by P u Il u tp u . In particular, in this projected subsystem the states useful 
for performing quantum computing are the reduced projected states ffg 00 (t). These can be 
obtained by using a projection operator Tr R P u Il u in Liouville space to act on the density 
operator p of the total system. Although p p s r ° 3 (t) is not a reduced density operator p$ (t), 
one can consider p p s ro:) (t) to be a generalization of ps (t) (i.e., when 11^ — ► P„, J2u Pv^v — ► 1, 
pp™? ^ _ TVb J2i> Pv^-vP (t) = Tr B p(t) = p s (t)). We argue that the information encoded 
in this projected subspace can be measured because the projected bases are orthogonal 
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and distinguishable in Liouville space - namely (j?J oj,s ^' oj ' s ^j = 0, for v ^ fi, because of 
Pfj,P u = in Liouville space too. 

(2) The phase shift induced by the eigenvalues may cost a unitary type error for quan- 
tum computation. Although this sort of error may be eliminated by developing standard 
quantum error correction schemes (such as multiqubit code' 24 '), the error recovery is not 
easy since phase error induced by the systematic phase shift may occur, at the same time, 
to many different clusters (inducing a phase or bit-flip errors in the encoded data). The 
increase in the number of phase errors for different qubits will cause a fast increase in 
the number of clusters, which likely leads to impractical implementations. On the other 
hand, the property of invariant eigenvectors in subspace provides the possibility of elimi- 
nating the unitary error by adjusting the appropriate time scale for the evolution operator 
to remain invariant under certain conditions; e.g., by choosing the time delay At to allow 
£ l ,e- < J'o JsS(T)dT P I , = j2 u e- l ti +At ( E » {T)+AE » {r) ) dT P u . Below, a concrete application of these 
concepts is presented. 



IV. COMPUTING SYSTEM 

We consider a two qubit quantum computing system S, consisting of spins Si and S 2 , 
such as those corresponding to two electrons around two 31 P confined in a germanium/silicon 
heterostructures in an electron spin- resonance transistor' 4 ' , or for two electrons confined in 
two quantum dots' 2,3 '. Ignoring the influence of the environment, the Hamiltonian may be 
written using the Heisenberg model as 

H s (t) = J(t)S 1 -S 2 , (13) 

where J (t) is the time-dependent exchange coupling parameter determined by the specific 
model considerations. In the case of spins of the two electrons (e.g., confined in two verti- 
cally (laterally) coupled quantum dots), J is the difference in the energies of two-electrons 
ground state, a spin singlet at zero magnetic field, and the lowest spin-triplet state; J is 
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also a function of the electric and magnetic field and the interdot distance! 25 ' 26 '. Using the 
relationship between Si ■ S 2 and the square of the sum of Si and S 2 , the eigenvalues and 
eigenvectors of Si • S 2 can be found from Si • S 2 = \ (S 2 — |) by 

|0!> = |11> , 

El s = \ |0 2 ) = |OO), (14) 

^> = ^(|01> + |10», 
E 2 s = -\ |04> = ^(|O1>-|1O». 

A quantum XOR gate is given by the sequence of operations' 2 ', Uxor — 
e t(n/2)sf e -i(iT/2)s^jji/2 e tnSfjji/2^ w h ere jj gw j s an (ideal) swap operator and determined gen- 
erally by an evolution operator U s (J (0) r) by adjusting the coupling time between the two 
spins in the evolution of the system. For the particular spin-spin coupling duration, r s where 
Jq s J (t) dt = 7r (mod 2n), U sw = U s (ir), the swap operator, is given by 

U sw = e-'Io 3 H ^ T)dT = £ e-* \4> 3 ) + \<f> A ) (0 4 | (15) 

and can exchange the quantum states of qubit 1 and 2. 

If we consider the influence of the environment, the non-ideal action of the swap operation 
must be considered because the influence of environment introduces decoherence. To treat 
this decoherence, one needs to understand the behavior of the evolution of the system. 
Here it is assumed that the environment consists of a set of harmonic-oscillators whose 
Hamiltonian is given by H B = J^k^kK^k, and the Hamiltonian coupling to the two qubit 
spin system is H int = XJ2k ( a l + a 1) {dk^t + 9k^k) , where b£, bk are bosonic operators 
for the &rth field mode, characterized by a generally complex coupling parameter which 
characterizes the case as being one of either independent or collective decoherence' 27 ' . The 
Hamiltonian operator for the total system is given by H (t) = H s (t) + H B + XH int , and the 
corresponding Schrodinger equation is ij^ip (t) = (Hs (t) + Hb + H int ) ip (t). Choosing the 
time-independent eigenprojectors of Hs (t) + H B as P u and Q u with Q v + P v — 1, with 

P v = \ip u ) (ip v \ = \ni---n k ---) \<f>j) (4>j\ (•••n k ---n 1 \ (16) 
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for v — (j, n\ • • • nk • • • ) and j — 1, • • • , 4; nk — 1,- • • , respectively, then the eigenprojectors 
U u (t) for the total Hamiltonian H (t) can be written in terms of the Heisenberg equation 
and satisfy the usual properties of projection operators. 



V. CODE CORRECTION IN SUBSPACE 

Using above subdynamics formulation, it is apparent that before the interaction occurs, 
the spectral decomposition for (Hg (£) + H B ) is the same as that for the total intermediate 
operator (t), i.e., 



H s (t) + Hb = (*) (17) 

4 



E E Ej,ni-n k - (*) Pj; 
j=l n\..n k — 



where i?° ni ... nj , ... (t) = (j — <5j4§) J (t) +Y,k ^fe^fe- But once interaction with the environment 
occurs, the eigenvalues change from £^ ni „, (i) to Ej ni ... nk ... (t) + AEj jnv .. nk ... (t), where we 
have made the following definitions: 

A-Ey, ni ... nfc ... (t) = (ni • • -7lfc- • • | Iln.iCj.,, . ..„,... (t) \4>j) |m • • -7lfc • • •) , (18) 

while the corresponding the eigenvectors of (t), | </>.,•) |ni • • • rik • • ■ ) remain invariant. The 
spectral decomposition of the intermediate operator (t) is thus different from H (t), 

©(*) = E E (n ■ a^, bj ,...(/!) /',„:...,„.... (19) 

j=l m-n k - 

This shows that there is no decoherence in the projected subspace for the stationary states 
of the system since the eigenvectors remain invariant, but for the evolutionary states, the 
change of the eigenvalues can introduce a type of unitary like error in the system evolution, 

U proj (t)) (20) 



4 

E E e 

j=i m-nk— 



- i /o( B °,ni-n fc -( T ) +A ^.»l-»*-( T )) 



dr 



ib proj (0) 
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As mentioned previously, the quantum XOR operator can be constructed from a se- 
quence of operations related to the ideal swap operator, adjusted by controlling the cou- 
pling time for the interaction between the two qubits without any influence from the envi- 
ronment. But if the effects of the environment are now included, the ideal swap operator 
changes to the non-ideal swap operator, owing to the unitary error. To cancel this type 
of decoherence, in terms of the subdynamics formulation, it is proposed that one allow 
the quantum logical operators to work on the projected subspaces. For example, if the 
quantum XOR operator previously introduced is considered, the ideal swap operator should 
be adjusted by controlling the coupling time between the two spins without considering 
interactions with the environment; then the ideal swap operator is given in a projected 
subspace as U sw (r s ) = e -^ s i- s ^B+is l sH B T 3 ) ^ w h ere a specific coupling duration t s is given 
by Jq 3 J (r) dr = n (mod27r). The unitary error is related to the non-ideal action of the 

swap operator which can be adjusted by the evolution operator in a projected subspace as 
U' sw (r s + At) = e -il? +At dr[J(r)s 1 -s 2 ^i B +i s m B +m nt c^ 1 ... nk ...(r)]^ The gpectral decomposition 

of U' sw (t s + At) can be expressed by adjusting the interaction time r s to r s + At, allowing 
the non-ideal swap operator to be equal to the ideal swap operator, 



n\—n k - 



U' sw (r s + At) = U sw 

m-n k - [j=l 

+ e ~ i lo^^ ' dT[-y{T)+Y, k ^ k n k +AE 3tni ... nk ...{t)} p 1 
n\—n k — \j=l 

This induces the integral equation for determining At which depends on j, n\ • • • 



(21) 



/ 

Jo 



r>+*t [J(t) f 3J(r) Arn , A f,forj = l,2,3, 

dr I -7^ " tv—T 1 + AE J,m-n k - (r)) = { (22) 

-f , for j = 4. 



For instance, if the exchange interaction coupling J (B, E,d) is a time-independent function 
of the external magnetic field B, electric field E and the interdot distance d as in the case 
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of vertically or laterally tunnel-coupled quantum dots, then At can be solved from Eq.(p2| 

as: 



At 



J(B,E,rf) _|_ j -3J(B,E,cQ 
4A_Ej > „ 1 ... nfe ... 4AB4 jni ... nfe .. 



+ r 



(23) 



where the concrete formula for J (B, E,cf), in the case of vertically or laterally tunnel-coupled 
quantum dots, can be found in refs.[25, 26]. Assuming, in this case, that the energy of the 
total system is uniformly distributed, i.e., 
J(B,E,d) -3J(B,E,d) 



4AE, 



4AE, 



const, for j = 1, 2, 3; n\, • • • , = 1, 2, 



4,m— rife- 



The shift energy can be obtained from Eq. fllS] 



ni---n fe ---) 



(24) 



(25) 



j'=l n' r --n' fc 



and the matrix of the creation operator is given by 



Cj',n' x —n' k — ;j,ni---nk--- 0Pj',n' 1 ---n' k --- ^j,nvrih— \ L Pj,ni---n k ---} > 



(26) 



which can be calculated from formula (|6|); for example, the first order C]^..,^ is divided 



by degenerate part and non-degenerate part: 



J',ni— n fe - 



(27) 



E 

^•••"fc- 

+ 5j' 4 



AE 



1 

j,ni-n k - 3',n' r -n' k 
2T EnJ-^- {Hint)j, n ^ 



Hint) 4' „' 



y,n,—n h ■■■ ;3,nvn k - 



Pj',n[-n' k - ;j,ni—n k - 



3,n v --n k - 



■n k - ;j,nv-n k - 



Pj,n' 1 ---n' k --- ;j,n\—n k — i for J 1 1, 2, 



0, for j = 3,4, 



which results in the second order AE^ ni ... nk ... by Eq.(|2 



AE 



int/j,n',—n', — ;j,n-L—n k — I ' 



j,m---n k - 



0, for j = 3,4, 



(28) 
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with the matrix elements given by 

(m • - • n k ■ 



(Hint) jini ... nk ... ; j> in ' l ... n ' k . 



(29) 



(<f>i\ H int \4> r ) 

A (Si d i - 5 2 ,f) E (5n k ,n> k +l9kyfn' k + 1 + Sn h ,n' k -igkifiE) > 

0, for j = 3,4, 

Eq.(P3"D shows that although the interaction introduces the sort of phase shift in the swap 
operator, this sort of phase shift can be cancelled by adjusting the coupling time between the 
two spins under the assumptions of homogeneous distribution of energy (i.e., owing to the 
invariance of the eigenvectors in the projected subspace). In the same way, the second-order 
projected states in the projected subspace are given by 

^j,ni-n k - fa) = {Pj,ni—rik— + ^j,ni—n k — — Dj,n 1 --n k --Cj t n 1 ... nk ...) {t)) 



(30) 



P 



(Sij + 5; 

,[2] 



2j, 



E {Hint) j, ni -n k - \j,n' 1 -n' k - Pj,ni—n k - \J,n{—n' k 



(S^ + 8 2 j) 



E E 



H, 



(AF [2] 

\ L ^ J - l 0,nx—nk — J "i "/. "l 

where the first order destruction operator is 



[H. 



j,m—nk—;j,n' 1 —n' k — \ lnt Jj,n"—n' k —;j,ni—n k —J 3,ni—nj, 



|V>(*)> 



j,ni-n k - 



(31) 



\-Hint)j >ni ... nk ...-j ;n ' l .,. n ' ... Pj^—n'k— ;j,m-n k -) f° r J 1)2, 



0, for j = 3,4, 

By above cancelling procedure, the evolution formula for second-order reduced projected 
density operator, before or after the interaction, remains the same and is given by 

$,nv»n h ----,j',nv--n k --- (*) ( 32 ) 

[2] 

H j,n 1 —n k —;j',n 1 —n k — V w / ' 

represent probability of existence for the 



where p J>r 



state 



(°)) (^n^-n^- (°) > notin g that f%ni°n k ---;j>,nv--n k --- (°) is S iven b Y initial de " 

.n,... ;/ ,„,..„ fc ... (0) = •/ ,•„■ (/';,,,..„,.../;>• (0) ® p B (0) /',■,/...„; ...)• The 
12 



w-oj[2] 



proj,S 



proj,S 



coupling condition, p^ ni J 



fidelity, which measures the decoherence between the initial state ip"P ro ^ (i) and the evolu- 
tion reduced project density operator ff ro ^' s (£) in the subspaces, is given by 



V j,j',ni—n k - 



(33) 




Z-^ rj,m— n k — ;i',ni— n*,— V w 7 



3,3 ,m---n k - 



j',nv-n k - l 



Pj,m—n k — Pj',ni-n k - 



3,ni-n k - 



3',ni-n k - 



which shows that there is no decoherence after the cancellation process. 



VI. THE DF PROJECTED SUBSPACE 

A completely decoherence-free condition (i.e., including no phase shift) in the projected 
subspace can also be determined, due to the definition of the intermediate operator of Eq.([|). 
That is 

H int C(t)=0. (34) 

In this case, eigenvalues for the intermediate operators in the projected subspaces, are the 
same as those for the total Hamiltonian in the total space, while the eigenvectors are the 
same as those of the unperturbed Hamiltonian. The projected subspaces spanned by the 
set of the eigenvectors of the intermediate operator are invariant for the evolutionary and 
stationary states of the system, and are closed for the intermediate operators. This allows 
quantum computing with coherent states such that no information leaks out from the pro- 
jected subspaces. Hence one can use the full power of DF subspaces to preserve coherence 
for full scale quantum computing. The key question then is how to ensure that Eq. ( |34| ) 
holds. In the Born-Markov approximation ( e.g., Lindblad equation approach) and the as- 
sumptions of symmetric or collective decoherence' 6,7 ' 13,24 ], if P M are partial eigenprojectors of 
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Hi n t (i-e., HintPfj, = EaPfj, , V degenerate a G set of complex number c) then it means having 
PyH int P^C v (t) P u = 0, for any v. But in the general case, without restrictions on the type 
of decoherence (e.g., a non-Markovian process and non-symmetric and non-collective deco- 
herence) how can one find the conditions for constructing the DF subspace? For generality, 
we propose a procedure to find the condition for the DF subspace starting directly from 
Eq.(|34|) by considering weak coupling between the system and the environment, without 
making any of the above assumptions. In fact, from Eq.(|4|) we have: 

Pj,ni---ni s ---Hi n tCj rll ... rlk ... (t) (35) 

4 1 
= X/ ~~rp~d 7+\ fpO (+\ T A TP 7^Pj,nv»n k -"HintPj',n' 1 --n' k -HintPj,m-n k - 

/ ln' r -n' k - ^3,nvn k - W ^j',n{-n/ k - \ l ) + ^-^j,ni-n k - IN 

A 2 

(dkgl+l ( n k + 1) + 9k-l9k n k) SjlPhm-nk- 



El„... (t) - I-y. „ (t) + AE 1>ni ... nh ... (t) , 

A 2 



, .i , , , ;/u f j\ T a TP ITS [9k9k+l ( n k + 1) + 9k-l9k n k) <hjl , i. ! .. -n i - 

= 0. 

This gives a DF condition which is a restriction on the bath operator by 

(9k9* k+ i (n k + 1) + gk-i9* k n k ) = 0. (36) 



(37) 



In one spectial case it is 

nk + 1 _ 9k-i9k 
nk 9kg* k+ i ' 

Under the DF conditions, e.g. with Eqs.(^) or fl3~7]) holding, Eq.(P reduces to 

d\(D proj (t)) ■ \ 

{Z "-=H {t) ^{t)) (38) 



dt 

with evolution of the projected state described by 

^proj ^ = f e -itH (t) \r>roj ^ (39) 

= f e -uH (t) ^2 (p v + Dv (o))- 1 ( Pv + Dv (0)) |p (0)) , 

V 

which enables the constructed projected subspaces to become a DF subspace spanned by 
the set of {Pj,n v ~n k -}i although the total space is subject to decoherence. The projected 
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subspace is closed with respect to the intermediate operator, 0„ (t) $ proj c $ pro - ? , for any 
projected state in the projected subspace $ pro - ? , (P„n„ (£) "0 (£)). Furthermore, if one general- 
izes the initial test projector P^ to be up (or down) triangular from the partial eigenprojector 
of H int , i-e. PuHintPfj, 7^ for v < /i, otherwise it is zero, then Eq.(|2]) gives C u (t) P v = 
which enables Eq.(]54]) to generally hold in the projected subspace without any approxi- 
mation. In this sense, if one can construct a triangular basis for Hi nt to span a projected 
subspace, then this space is DF although the total space is decoherent. 

VII. REALIZATION OF A DF SUBSPACE USING TRIANGULATION 

As a starting point we consider a necessary and sufficient condition for DF behavior in 
the projected subspace. This can be determined from the fact that the interaction part of 
the intermediate operator 9 (t) is zero or that the evolution of the projected (reduced) state 
is independent of the interaction part of the Liouvillian. That is 

L nt + (Lf (t) + Lf + L\ nt ) C (t) = 0. (40) 

If the constructed projectors P v and Q u are triangular with respect to Lf (t) + 
Lf + Lf* (i.e., P u L mt P^ ^ o, P v (Lf (t) + Lf ) P„ ^ for v < fi (up- 
triangular) or v > jj (down-triangular)), then L nt + (Lf (t) + Lf + Lf 1 *) C (t) = L nt + 
Yn, P» (Lf (t) + Lf + Lf*) Qu Ev(t) _Q vmQv Qu (Lf (t) + Lf + L\ nt ) P v = 0. This enables 
the evolution of the reduced projected density operator to be independent of the interaction 
part of the total Liouvillian in the projected Liouville subspace, 

(t) = £ Ttb [{Pv + Dv {t)) p {t)] = ^ Tr B (e~< dt ' lf0{t,) P u p (t )) , (41) 

where If (t) is a uth eigenvalue of 

Go (t) in Liouvillian representation. The corresponding mixed-state fidelity is equal to 1 
(i.e., F(t) = Tr s {ff Toj ' s (0)EvTr B (e~*^ dt ' l "° {t,) P„p (t )\ J = 1). This shows that there 
is no decoherence introduced by the environment in the projected subspace, although the 
total system is subject to decoherence introduced by the environment. The formulation for 
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this DF projected subspace is exact (i.e., there are no approximations including the Born- 
Markov approximation). The method is general. Thus for any combined system, one may 
construct a triangular basis for the total Hamiltonian to span a projected subspace in which 
the evolution of the projected density operator is independent of the interaction part of the 
total Hamiltonian, with a fidelity of 1. Furthermore, one can construct a partial triangular 
basis for the computing system or bath, by choosing a partial diagonal basis for the left 
part of total system to form a DF projected subspace, based on the subdynamical Liouville 
equation. Thus triangulation procedure places no restriction on the type of coupling between 
the system and the environment (i.e., Markovian or collective decoherence) . 



VIII. AN EXAMPLE 

To illustrate the above method, we revisit the example of the two qubit quantum com- 
puting system S, consisting of spins Si and S 2 . For example, this includes the two elec- 
trons around two 31 P confined in a germanium/silicon heterostructures in an electron spin- 
resonance transistor^ ; or two electrons confined in two quantum dots' 2 '. Ignoring the in- 
fluence of environment, the Hamiltonian can be written using the Heisenberg model as 
H s (t) = J (t) Si • S 2 , where J (t) is the time-dependent exchange coupling parameter deter- 
mined by the specific model considerations. The coupling to the environment is assumed to 
be described by a Caldeira-Leggett-type model, consisting of a set of harmonic oscillators 
coupled linearly to S by H int = A J2l=i °~k • bfc, where V k = g% (a a>kj + a+ fej ) is a fluctua- 
tion quantum field. The Hamiltonian of bath is given by H B = J2 a ' a a,kj a a,kj , where a+ fe ■ 
( a a,kj) are bosonic creation (annihilation) operators with j = x,y, z and lu%> are the cor- 
responding frequencies with spectral distribution function J k j (00) = n J2 a (fl'a' 7 ) ^ ( u ~ w a) 

[28] _ 

We show how to triangulate the partial basis of the bath since, in principle, triangulating 
the partial basis of the system follows the same approach. The key point is that, in general, 
the triangular property of L results in P V LQ V C V P V being in the SKE. 
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The matrix for Hb + H int , with respect to a cuth element of the basis of B, in a repeat 



subspace, ( 



n 



a,k / ' 



<fc + l)> (<fc|, (<fc + 1 |)> 



can be written as 



M 



^'<* 9a J \/K,k + 



V 



,/. - ' "n't"/,./, • 1 



Using a similarity transformation to up-triangulate M, 



( 



M 



tri 



ad — be 
1 



d -b 
— c a 
( 



( 



a lb a,k 



kjj n 3 



a,k T 1 W a (<k + 1 



\ c 7 



ad — be 



adw k J^ k (d 2 - 5 2 ) ^V<* + 1 - 
^ (< fe + 2) + adu k J + 1 



we have 



- bcug « fc + l) + (dc - oft) ^V<* + 1 = °> 
(-c 2 + a 2 ) <?*V< fe + 1 + aav k J = 0, 



which gives 



± 



a;, 



+ 1 c = 7 Q)fci c, 



kj 
9<* 



-b = Ca,kjb- 



(42) 



(43) 



(44) 



This determines the triangular basis for + if int in the repeat subspace as 
( <j)n <8> { ra a,t}) , ® 0n )• For example, by choosing c = 1 and 6 = — 1 we have 

\ 



M 



tri 



-J*,kj(a,kjiOa j K,k (C,kj - l) 9a\/ni,k + 1 ~ Ca.krf 



(<* + 2) - -y a ,kjta,kjW$ (< fc + 1) 



(45) 



defining, 



a, A; 



a, A; 



7a, fc? 



n 



a,k 



ni.k + 1 



ct,fc 



Ca,kj 



nik + 1 



(46) 
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and 



rv 



j,+ 

a,k 

j- 
a,k 



1 



~'~fa,kjCa,kj 1 



-c 



a,kj 



7V 



— ( n 



L a± 



\ 



n 



a.k 



+ 7a,fcj (<fe + 1 



(47) 



with jn^j = f^i'fc • ■ ' ' ')• Then we can construct a triangular basis for the total 
Liouvillian L (t) by {||0 n ® {ng}) ({ng} ® 0m|), (|0m ® {njj}) ({<fc} ® 0n||}- This 
gives triangulation of matrix of Lb + Pint- By means of the new triangular basis one can 
define the projectors as P u = 4> n ® {n 3 ^}) ({^'|) ® 0m ) ( 0m <8> {^'l}) ® 0n 

and (3„ = 1 — P v . Taking into account the definitions of the creation and destruction opera- 
tors, one has P v (Pf + Lf *) P^C,, (t) P, ~ P, (Pf + Lf *) P M (Lf + Lf *) P, — ► for > v 
(up-triangular property). Therefore one finds the DF condition: Lq™' + (Lf + Lf *) C (t) = 
with (t) = O (t) = Lq + Lq. Under this DF condition, the SKE reduces to 



.d\p pro Ht)) 



dt 

= O (t) ff r ° 3 ' (t) . The evolution of the reduced projected state, in the initial decoupling 
condition, becomes to Eq.fl41|). This shows that p p J OJ (t) is independent on L int . Therefore 
there is no decoherence introduced by L int in the projected subsystem although the total 
system is subject to the decoherence introduced by L int . In this projected subsystem, the 
quantum Control-Not logic operation is still given by a sequence of operations and the swap 
operator U sw remains invariant before and after the interaction, and is given by 



Tr B fe- < ^'* e(T) *'p proj '(0) 



Tr B fe~ l r^W^y^o) 

Tr B (fe~^o S H o(r)dr pproj ^ fj ft' H§(t)<1t 

Tr B (U swf f ^ (0) P s -J) , 



(48) 



where 



Ps 



EE' 

a \n=l 



W (la 



+ e l 



1«>(1 £ 



(49) 



The following table clarifies the difference between the total space and the projected 
subspace constructed by the triangulation: 
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Initial State 


Evolutionary State 


Fidelity 


Subspace 


Tr B PvPs (0) ® p B (0) 


Tr B fe- l S L ^ dt P v p s (0) ® p B (0) 


f = l 


Total Space 


Tr BPs (0) <g> p B (0) 


Tr s T e" l /( L o W+if +^f*(*))* ps (o) ® Ps (0) 


0</<l 



If one assumes that in the Schrodinger picture, an initial state of 5 is 
05(0) = |(|TI) + |IT)) - \ (ITT) + III)), then by choosing a projector P v as 
\ (IT!) + lit)) ® { n «,fc}) ({ n a,fc} ® | (IT!) + IIT)) , after taking into account the obtained 
triangular basis, an initial state in the DF P^-projected subspace is given by TtbP v <1>s (0) = 
HIT!) + IIT)) under the initial decoupling condition. This correspondence is realized by the 
projection: <ps (0) — ► P„4>s (0). The advantage of using P^(f>s (0) compared with (f>s (0), is 
that the evolution of P^<f>s (0) is independent of the interaction between S and B in the pro- 
jected subspace. Therefore, one can encode information onto P„4>s (0) to perform quantum 
computing (or quantum communication), protecting against decoherence introduced by the 
interaction between S and B. 

Finally, it should be noted that it is not necessary that the evolution operator Uq in 
the projected subspace is unitary, because a quantum computing system projected in a DF 
subspace may be an open quantum system, obeying the semigroup evolution rules^ 17 !. In 
this open system, self-adjoint operators and unitary evolution groups are not intrinsically 
necessary to govern quantum computation. Quantum computation can then be performed in 
a more general functional space, such as RHS, rather than just Hilbert space. The projected 
state ^ proj may exist in the test space $?' ro -?, which is a dense subspace of the Hilbert space 
q-(proo constructed by J2 U Pv^-v^i-, representing the physical states which can be prepared in an 
actual experiment. Its adjoint ^} pro ^ lies in the dual space (<3> proj ) x , representing a procedure 
that associates with each state a number, while preserving the linear structure which results 
from the superposition principle, i.e., the triplet structure $P r °.? c. 7i pro ^ c [<^ pro ^) x I 23 l. This 
is a RHS structure which facilitates describing irreversible processes like decoherence and 
dissipation due to interaction with the environment. In this space the evolution of the states 
are permitted to be time asymmetric, providing a framework for describing the irreversibility 
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of practical open systems. This irreversibility does not change quantum reversible logical 
operations to quantum irreversible logical operation in the quantum universal Controlled-Not 
logical gate. To appreciate this one must distinguish between irreversibility of a quantum 
logical operation, introduced by the structure of logical gate, and irreversibility of the process 
induced by interactions with the environment. Reversible computation means reversible 
logical operations on the structure of the logical gate. In this sense, quantum computing in 
RHS is compatible with reversible quantum logical operations and permits computing any 
reversible function, although irreversible processes do in fact exist. 

IX. CONCLUSION 

In conclusion, a subdynamics based formulation in the Schrodinger picture was presented 
for an open quantum system. Based on the subdynamical kinetic equation for an open quan- 
tum system, a proposal for quantum computing in the projected subspaces is developed. The 
eigenvectors of the intermediate operator in this subspace were shown to remain invariant 
before and after interaction with the environment, while the eigenvalues of the intermediate 
operator in this subspace change. The fidelity of mixed states in the projected subspace 
is 1 which means that the constructed projected subspace is definitely DF. This reveals a 
universal property for any system plus reservoir: one can construct a DF projected subspace 
by using eigenprojectors of the free part of the Hamiltonian in which the encoded states 
are projected states, which are themselves determined by relevant formulae of subdynamics. 
On the other hand, changes of the eigenvalues after interaction may introduce a type of 
unitary error in the ideal swap operator. This sort of error can be cancelled by adjusting the 
interaction coupling time between two spins in the subspace, since the eigenvectors remain 
invariant, although decoherence exists in the total space for the total system. Finally, the 
general case for completely DF behavior in the projected subspaces was discussed, and it 
was shown that using a general DF condition (i.e., the second term of the subdynamics 
kinetic equation is zero), one can find a condition to allow the constructed projected sub- 
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space to be DF. This reveals that this condition is a necessary and sufficient condition for 
constructing a DF projected subspace. We wish to emphasize two points here: (1) the con- 
structed projected subspace is spanned by a set of {ffj ^ (t)}, which is closed with respect 
to the intermediate operator. Indeed, (t) ff roj (t) = £„ P »L (t) (P v + C v (t)) Yi v (t) p (t) = 
J2u if (t) P^ oj (t), where If (t) is a vth eigenvalue of 6 (t) in Liouvillian representation. (2) 
The Born-Markov assumption and various other types of restrictions for DF behavior does 
not need to be made, indeed as we show a general approach can be used. 
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